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3 , Abstract. The recently formulated framework of anisotropic and dissipative 

S . hydrodynamics (ADHYDRO) is used to describe non-boost-invariant motion of 

the fluid created at the early stages of heavy-ion collisions. Very strong initial 
asymmetries of pressure are reduced by the entropy production processes. By the 
^^JH ' appropriate choice of the form of the entropy source we can describe isotropization 

t^ times of about 1 fm, which agrees with the common expectations that already at 

Cn , such times the perfect-fluid hydrodynamics may be applied. Our previous results 

^D i are generalized by including the realistic equation of state as the limit of the 

\0 ' isotropization processes. 
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1. Introduction 

In a recent paper we have introduced the new framework of highly- Anisotropic and 
strongly-Dissipative HYDROdynamics (ADHYDRO) [1 . This framework allows for 
studies of the systems which are locally highly anisotropic. Such situations happen 
at the early stages of relativistic heavy-ion collisions where, as suggested by many 
microscopic calculations [21 [3l H], the longitudinal pressure is much lower than the 
transverse pressure [j. 

In Ref. [T] our approach has been used to analyze the behavior of matter created 
in the central region of relativistic heavy-ion collisions, where the assumption of boost- 
invariance is acceptable. We have shown how our equations describe the isotropization 
of pressure and how the system approaches the regime governed by the perfect-fluid 
hydrodynamics . 

In this paper we relax the assumption of boost-invariance and consider a general 
one-dimensional expansion along the beam axis (at the early stages of collisions the 
longitudinal motion dominates). We derive equations describing non-boost-invariant 
motion of anisotropic fluid and show how the initial anisotropic behavior transforms 
smoothly into the locally isotropic expansion. 

In the second part of the paper we present a scheme implying that the 
isotropization transition leads to a desired equation of state in the perfect-fluid stage. 
In this way we construct an effective model which describes several stages of heavy-ion 
collisions in a very concise way: an early anisotropic phase (Glasma, color-flux tubes), 
transition to the isotropic phase with the QCD equation of state, perfect-fluid stage 
(sQGP), and the phase transition to the hadron gas. 

We note that similar physical problems connected with the non-boost-invariant 
motion of the fluid formed at the early stages of relativistic heavy-ion collisions were 
studied in Refs. 01^, see also [TUl [TTl [T^]. In Ref. j5j the effects of the equation of 
state and the initial conditions were studied, while in Ref. [6] the effects of the shear 
viscosity were taken additionally into account. In contrast to Refs. [5l|6] our approach 
concentrates on highly asymmetric stages of the evolution, see |71[8]. We also note that 
our formulation of ADHYDRO was followed by the paper by Martinez and Strickland 
[5], where a similar framework of anisotropic hydrodynamics was introduced. We 
discuss the connection between our work and Ref. ^9, in the Appendix. 

2. Anisotropic hydrodynamics 

In this Section we recapitulate the main assumptions of ADHYDRO, that have been 
recently introduced in [J. Our starting point is the following form of the energy- 
momentum tensor, 

T^'' = (e + p^) U^'U" - P± g^"" - (P_l - PwW^'V. (1) 

Here e, P±, and Pii are the energy density, transverse pressure, and longitudinal 
pressure, respectively. For isotropic fluid, where the two pressures are equal, 
Pj_ = P|| = P, we recover the form of the energy-momentum tensor of the perfect- 
fluid hydrodynamics. The four-vector If^ in ([T]) represents the hydrodynamic flow, 

C/'^ = 7(l,«.,«y,^'.), 7 = (l-«')^'^', (2) 

I The directions are defined with respect to the beam axis 



while y defines the direction of the longitudinal axis that plays a special role due to 
the initial geometry of the collision, 

^^=7.K, 0,0,1), ^. = {\-vl)-^l\ (3) 

The four-vectors U^ and V^ satisfy the following normalization conditions: 

C/2 = 1, T/2 ^ -1, [/-F^O. (4) 

In the local-rest-frame (LRF) of the fluid element we have U^ = (1,0,0,0) and 
V^ = (0, 0, 0, 1), and the energy-momentum tensor has a simple diagonal structure |§|, 
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Hence, as expected, the formula ([T]) allows for different pressures in the longitudinal 
and transverse directions. Anisotropics of the form ([5]) are present in the case where 
the particles are formed by the decays of color strings and in the theory of Color Glass 
Condensate [H [H [H] . 

Besides the energy- momentum tensor ([1]), we introduce the entropy flux 

o-^ = crC/^, (6) 

where a is the entropy density. We assume that e and a are functions of Pi and P\\ . 
In particular, for massless partons the condition T^^ — gives 

£ = 2P^-|-P„. (7) 

The form of the energy-momentum tensor ([1]) resembles the form used in 
relativistic magnetohydrodynamics [TBI El ■ In that case the anisotropy is induced 
by the presence of the magnetic field. At the early stages of heavy-ion collisions we 
have similar situation — there exist strong color magnetic and electric longitudinal 
fields (Glasma following CGC [H]) which polarize the medium. In this paper, however, 
we do not include explicitly the effects of the fields. Such analysis is left for a separate 
study. 

The space-time evolution of the system is governed by the equations expressing 
the energy-momentum conservation and the entropy growth, 

a^fi^ -S. (9) 

The function S represents the entropy source. The form of E must be treated as an 
assumption that defines the dynamics of the anisotropic fluid. It is natural to assume 
that E > and S = for Pj_ = Py . In this way, in the case where the two pressures 
are equal, the structure of the perfect-fluid hydrodynamics is recovered. 

We treat E as a function of Pi and P|| . In this way, Eqs. ([5]) and © form a 
closed system of 5 equations for 5 unknown functions: three components of the fluid 
velocity, Pj_, and P||. The projections of Eq. ^ on Uy and Vy yield 

§ The Lorentz transformation which leads from the center-of-mass frame, where Eqs. l[2j and ((Sjl are 
vaUd, to LRF, where we have V^ = (1, 0, 0, 0) and V^ = (0, 0, 0, 1), consists of a Lorentz boost along 
the z-axis, a rotation around the z-axis, and a boost in the direction of the remaining transverse part 
of V. In this way the initial longitudinal direction defined by the spatial part of V^ is left unchanged. 






(e + P^) d^U'' + {Pi_ - P||) U^V^d^V, 
{P\\ - P±) d^V^ + (e + P±) V.U^d^^U". 



(10) 

(11) 



3. Anisotropy parameter x 

Our earlier studies show [Ij [T9j [20] that it is useful to switch from P± and P\\ to the 
two new variables: the entropy density a and the anisotropy parameter x. In this case 
we may write 



(12) 
P± = I ^-^ 1 ^-^^^+xR'(x) , (13) 

Pii = I ^^^ 1 ^^^ - 2xR'(x) . (14) 
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The function R{x) is defined by the formula [2(1 1 
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(15) 



and go is the degeneracy factor connected with internal quantum numbers of particles 
that form the fluid. In our calculations we take go = 16 — the initial system is formed 
most likely from gluons. The symbol R'{x) denotes the derivative of R{x) with respect 
to a:, for X = 1 we have R'{1) = and the two pressures are equal. 

The physical interpretation of x and details describing manipulations that lead to 
(fT2l) - (fT5|) are given in [U |20]. Here we only note that the ratio P±/P\\ is a monotonic 
function of x with Pi > P|| if x > 1 (to a good approximation P^/P^^ = a;"^/^, and 
P± =P|| if 2; = 1). 

4. One-dimensional non-boost-invariant motion 

For a purely one dimensional motion we may introduce the following parameteriza- 
tions: 



W = {coshd{T,r]),0,0,smhd{T,r])) 



and 



V = {smhd{T,ri),O,O,cosh'0{T,r])) , 
where d{T,ri) is the fluid rapidity depending on the (longitudinal) proper time 



= Vt^ 



(16) 
(17) 
(18) 



II Note that for x < 1 the function (arctan V^ — ^)/V^ ^ 1 should be replaced by 
(arctanhx/l — x)/V^ ~ ^ 



and the space-time rapidity 

1, t + z 
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(19) 



Equations ([16]) and p7|) satisfy automatically the normalization conditions (|4]). We 
also have 

d sinh(^ — rj) d 



U^'d^ = coiih{-d-r])— + 



(20) 



dr T drj 

The expressions similar to (l20l) may be derived for d^JJ^, V^df^, and df_iV^. Using 
them in Eqs. ((TU]) and ([TT|) gives 



d tanh(i9 — rj) d 
dr T 



tanh('i9 — rj) 
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dr rdrj 
Similarly, from Eq. ^ one obtains 
d tanh('i9 — rj) d 
dr T drj 

The structure of Eqs. (PT|) ~(P5 |) is very much similar to the equations studied in Refs. 
0. 
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d tanh(^ -rj) d' 
dr T drj 
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cosh(z9 — rj) 



5. Entropy source 

In the boost-invariant case the functions e, Pu, P±, and a do not depend on rj, and 
d — rj. In this case Eq. ([2^ is automatically fulfilled, while Eqs. (PT|) and (P^ are 
reduced to 



de 
lb- 
da 
dl 
Equations dMl and 
source E was used 



-'-^±^. (24) 

T 

- = E. (25) 

r 

were studied in [T] where the following form of the entropy 
(1 - V^)2 a 



(26) 
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Here Tcq is a timescale parameter. The form (pS)) guarantees that E > and 
E(cr, a; = 1) = 0. More arguments for such a particular form of E are given in [1]. 
They connect x with the ratio of the transverse and longitudinal temperatures. 

It may be shown that for small deviations from the equilibrium Eq. (|26p leads 
to quadratic dependence of the entropy source on the variable ^ = 1 — a; and 
this dependence is compatible with the Israel-Stewart theory (where the entropy 
production depends on the viscous stress squared) and with the Martinez-Strickland 
model [9]. 

We have to stress, however, that the structure of the entropy source is an external 
input for the anisotropic hydrodynamics. It would be interesting to obtain any hints 
about E for large x from the microscopic models of particle production or from the 
AdS/CFT correspondence. 



6. Isotropization 



If £ and P|| are expressed in terms of a and a;, Eqs. ([2T|) -(|23 |) become three equations 
for three unknown functions: <j, x, and d: 
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(28) 
(29) 



Here we have introduced the ratio of the longitudinal pressure and energy density, 
which equals 1/3 for the isotropic system with a; = 1, 



h{x) 



a, m-\. 

Substituting ([^ into (P7|) and using the relation 



we find 



l + h{x) = 
id 
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Here we use the notation introduced in [T], 

R{x) (1 - Vi)' 



H{x) 



R'{x) 



(33) 



In the case where d{T,rj) sa rj^ Eq. (I32p may be approximated by the ordinary 
differential equation derived in [T] 



dx 2x AH{x) 
dr T 3T( 



(34) 
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Thus, in this case we may immediately use our previous results to argue that a: — > 1 for 
T >• "Tcq- If the condition ^{T^rj) « 77 is not fulfilled, one should study Eqs. (P7)) -(P ^ 
numerically. 



7. Initial conditions 

The initial conditions for a non-boost-invariant evolution are defined by three 
functions: (t(to,?7), i?(to,?7), and x{To,r]), where tq is the initial proper time. In 
the numerical calculations we assume tq = 0.2 fm. The initial entropy profile is taken 
in the form [HJ [22] 



a{To,r]) = croexp 



9{\r,\~A^) 



(35) 



2((5r;)2 

where 9 is the step function, the parameter A77 defines the half width of the initial 
plateau in spacetime rapidity, and Si] defines the half width of the Gaussian tails on 
both sides of the plateau. By the appropriate changes of the parameters Arj and Srj 
we may vary between the boost-invariant-like and Gaussian-like initial conditions. To 
match the rapidity distribution measured by BRAHMS ^23^ we use the values: A77 = 1 
and Sr] = 1.3. The value of the initial central entropy density is obtained from the 
condition that the initial energy density is 100 GcV/fm'^, 

/ 2 \ 4/3 

£0 = 100 GeV/fm3 = ( ^ I R{xo). (36) 



Here, xq is the initial value of the anisotropy parameter at 77 = 0. In the present 
calculations we assume that the initial profile of x is constant, 

a;(To,r/)=xo = 100. (37) 

This value of xq is somewhat arbitrary but it expresses the opinion that the initial 
longitudinal pressure is much smaller than the initial transverse pressure. 

The value of the timescale parameter Toq has been set equal to 0.25 fm. Finally, 
in agreement with other hydrodynamic calculations which address the problem of the 
longitudinal expansion, we take the initial fluid rapidity profile in the form 

Hto,v)^V- (38) 

Of course, for the boost-invariant motion we have always 'd = r]. In our case, due to 
the presence of the longitudinal gradients determined by the entropy profile (|35p , we 
expect I'd] > I77I. 

8. Results 

The results of our numerical calculations with the initial conditions specified in the 
previous Section are presented in Figs. [1] and [21 The upper part of Fig. [1] shows the 
time evolution of the anisotropy parameter x for two different values of the space-time 
rapidity: 77 = (solid line) and 77 = 6 (dashed line). We observe the initial rapid 
decrease of x (from 100 to about 5) followed by a longer evolution where x approaches 
unity. This behavior is very weakly dependent on rapidity and very much similar to 
that observed previously in the boost-invariant case. 

The lower part of Fig. [T] shows the time evolution of the ratio P^^/P±. The time 
dependence of this ratio reflects the time dependence of the anisotropy parameter 
X. Again, we observe a rapid change at the very beginning of the evolution that is 
followed by a more moderate approach to unity. For r > 2 fm the system is practically 
isotropic. Clearly, by changing the value of Teq we may speed up or slow down the 
process of isotropization. The most common perfect-fluid description of the data 
requires that the stage described by perfect hydrodynamics starts at about 1 fm or 




Figure 1. Upper part: Time evolution of ttie anisotropy parameter x for two 
different values of the space-time rapidity: t; = (solid line) and rj = 6 (dashed 
line). Lower part: time evolution of the ratio Pa/Pj_ for rj = (solid line) and 
r) = 6 (dashed line). 



even earlier. Therefore, we think that the time evolution obtained in our approach 
may be regarded as realistic if the systems becomes isotropic at about 1 fm. The 
results shown in Figs. [1] and [5] suggest that this condition may be achieved if r^q < 0.1 
fm. Such extremely short time scales indicate very strong interactions in the early 
anisotropic plasma. 

The upper part of Fig. [5] shows the difference between the fluid rapidity i) and the 
space-time rapidity 77 for t = tq (solid line), r = 1.2ro (dashed line), and t = 10 fm 
(dotted line). In the boost-invariant case t? = ry. In the present calculation, due to 
the extra pressure gradients determined by the entropy profile, the values of the fluid 
rapidity become larger than the space-time rapidity (for 77 > 0). 

The lower part of Fig. [2]shows the entropy density profiles in rapidity for different 
values of the proper time. The solid line shows the initial entropy profile psp with 
the central entropy density ao determined by the condition ((36)) . The thick dashed 
line shows the entropy density profile for the proper time t = 1.2to. Since the 
initial dynamics is dissipative, the entropy is produced in the early stage, and the 
thick dashed line is placed above the solid line. For comparison, the thin dashed 
line shows the entropy profile expected from the boost-invariant isentropic expansion 
where cr(r, rj) = a(To, vi)t^/t. The thin dashed line is placed below the solid line, since 
the entropy density decreases in this case due to the longitudinal expansion of the 



It 
I 



T= 1.2 To 



--■^""'''•""" 



T = 10 fm 




Figure 2. Upper part: difference between the fluid rapidity ■!? and the space-time 
rapidity t] for t = tq (solid line), r = 1.2to (dashed line), and t = 10 fm (dotted 
line). Lower part: profiles of the entropy density plotted for different values of 
the proper time. The results of our numerical calculations are represented by the 
thick lines. The results obtained from the simple scaling (t(t, 77) = o-(ro,»7)ro/r 
are represented by the thin lines. 



fluid. By comparing the maximal values of the thick and dashed lines we can make 
the estimate of the relative entropy production due to the dissipative effects present 
in our approach. From this simple calculation we find that the entropy increases by 
about 70%. 



9. Transition to quark-gluon plasma 

The concept of AD HYDRO is motivated by the microscopic picture of the anisotropic 
parton distribution function [U [5D]. As the system becomes isotropic, the parton 
distribution becomes also isotropic (locally, in the momentum space). This, however, 
does not automatically mean that the isotropic system is described by the realistic 
equation of state. The Boltzmann distribution used in Refs. (HHO] may be inadequate 
for description of the strongly interacting quark-gluon plasma. To connect the process 
of isotropization with the process of formation of the equilibrated quark-gluon plasma 
we may consider the following ansatz for the energy density and pressure 
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Figure 3. The same as Fig. [T]but with Eqs. |[T2)| - ||14| | replaced by Eqs. lf39)l - 
(ETJ. 



eia,x) ^'j^^Ri:,), 



<j go 

P±{(T,X) =Pqgp(cr) — 

90 

Pjl(cr,x) =Pqgp(cr) — 
50 



i?(a;) 

3 
R{x) 



xR {x 
2xR'{x) 



(39) 
(40) 
(41) 



Here, the functions £qgp(cr)/3 and Pqgp((T) describe the reahstic equation of state 
constructed in ^. They appear in Eqs. (I5^ -(PT |) in the very much analogous way 
as the term a^f^ in Eqs. ([T^ - p^ . Although eqgp(cr) and Pqgp(cr) were calculated 
for equilibrium, they are used by us to calculate the energy density and pressure for 
an arbitrary state characterized by the entropy a and the anisotropy parameter x. 
In order to take into account the effects connected with the anisotropy, we multiply 
eqgp(cr)/3 and Pqgp((7) by the same combinations of the function R{x) as in Eqs. (fT2|) 



p4p. The factor vr^/go has been introduced to guarantee that e{a,x = 1) = eqgp(cr), 
and similarly P\\ (cr, x ^ 1) = P^{a,x ~ 1) ~ Pqgp(o-). 



We should admit that Eqs. ((39)) - (l4T]) have no direct microscopic explanation. 
The main argument for using this form is that it has two attractive limits and may 
be used as an interpolating formula. At the very early stage, the system consists 
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Figure 4. 


The same as Fig. [2] but again w 



ith Eqs. |(T2)| - ||14|I replaced by Eqs. 



most likely of massless partons and Eqs. (|39| - pT|) approach Eqs. (|12p - (|T4)) . since the 
realistic equation of state approaches the Stefan-Boltzmann limit (although deviations 
from this limit are noticeable). On the other hand, if the system becomes isotropic, 
X — 7> 1, R'{x) — 7> 0, and Eqs. ([5^ -(PT |) are reduced to the equation of state used 
in standard hydrodynamics. From the statistical point of view, the initial system of 
partons is always treated classically and R{x) is given by (fTS)) . Inclusion of quantum 
statistics and consequences of other model assumptions will be discussed in a separate 
paper. 

In Figs. [3] and H] we show the results of the numerical calculations where Eqs. 
P9 P " (PT|) were substituted in Eqs. d^Tt - ipS]) . By this substitution we have obtained 
the system of three equations for ct, x, and d, which is quite analogous to Eqs. (j27p - 
P9p. The values of the initial parameters and the initial profiles have not been changed 
except for the initial value of the entropy density, which has been modified to keep 
the initial energy density unchanged, 

-" i?(xo) 



£o = 100GeV/fm3 = eqgp((To)- 



(42) 



50 3 

The overall space-time evolution of the system shown in Figs. [3] and |4] resembles 
very much the evolution shown earlier in Figs. [I]and[2l In particular, the time changes 
of X and the ratio of pressures are practically the same. The small differences are 
visible in the evolution of the entropy density (mainly due to a different normalization 



12 

caused by fixing the initial energy density rather than the entropy density) and the 
fluid rapidity. The difference -d — ij becomes now a more complicated function of 7], 
the effect noticed also in Ref. [B] and connected with the use of the equation of state 
from Ref. [H]. 

Clearly, an attractive feature of using Eqs. ([39)) - (|4T|) instead of Eqs. (fT2|) - ((T4| 
is that the system reaching equilibrium is described by the realistic QCD equation 
of state. In this way, the intermediate hydrodynamic evolution agrees with our 
expectations and no extra construction (like, e.g., the Landau matching condition 
used in |25| ) is necessary to change from one thermodynamic regime to another one. 
On the other hand, the connection of Eqs. (|39)) -(|4T |) to any microscopic picture is 
lacking, since the realistic QCD equation of state cannot be simply interpreted in a 
quasi-particle picture. Another attractive feature Eqs. (pn)) - (PT|) is that they include 
the phase transition from the quark-gluon plasma to the hadron gas. Thus, using Eqs. 
P9 p ~ (|4ip we may describe many stages of heavy-ion collisions in a single approach. 

10. Conclusions 

Our results confirm very attractive features of ADHYDRO as an effective model for 
early stages of heavy-ion collisions. In this paper we have shown that the tendency to 
approach isotropy is not changed by non-boost-invariant profiles. Very strong initial 
asymmetries of pressure are reduced by the entropy production processes. By the 
appropriate choice of the form of the entropy source and the value of the timescale 
parameter, we have succeeded to reproduce the realistic times for the isotropization 
expected in heavy-ion collisions. We have also generalized our previous result by 
including the realistic equation of state as the limit of the isotropization processes. 

From the practical point of view, our formalism based on the anisotropic 
distribution function allows for the determination of the space-time evolution of color- 
neutral anisotropic distributions, which may be used, for example, as background 
distributions in the analysis of the plasma instabilities. 

The results obtained within ADHYDRO in this paper, and also in our previous 
publication, open perspectives for further applications of the model. In particular, 
the model may be applied in the 2+1 case (longitudinal boost-invariance without 
cylindrical symmetry) or in the general 3-1-1 case. In the 2-1-1 case it may be treated 
as a generalization of the framework where transverse hydrodynamics is combined 
with the perfect- fluid hydrodynamics by the Landau matching conditions [3S]. In 
ADHYDRO the change from the system where only transverse degrees of freedom are 
thermalized to the full local equilibrium is described in the continuous way. 

11. Appendix 

Our formulation of ADHYDRO was followed by the paper by Martinez and Strickland 
[5] where similar ideas appeared. It is interesting to show that the two approaches 
agree for small deviations from equilibrium in the boost-invariant case. In the 
following, it is useful to use the variable 

^ = x-l. (43) 

Equation (20) from [9 connects ^ andphaid, 

1 d^ 2 6 dphai 



1 + ^ dr T phard dr 



2r 



l-7^3/4(e)vTT^ . (44) 
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Figure 5. Functions H{^) and (15/2)Hjvf5(5) used in the present paper and in 
Ref. [9]. For small ^ they are both well approximated by the linear form (45/16)^. 



Here, phard is related to the average momentum in the parton distribution function 
and TZ{£,) is defined by the formula 
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arctan y^ 
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(45) 



For 1^1 < 1 we expand 7?.(^) around zero. Keeping the leading terms in ^ one obtains 
from (iMl) 



1 d£, 



6 dpha 



re 

15 ' 



(46) 



l + ^dr T phard dr 
In order to connect this result with our approach, we relate the non-equilibrium 
entropy density a to phard [S] ■, 

a = Apl,,^x-^'\ (47) 

Here A is an irrelevant constant. Inserting Eq. (j47| into Eq. (|25|) and using our 
ansatz for E wc find 

Idx ^2 ^ 6 dphard _ 2 (1 ~ Vi)^ ^^ {X - If 

X dr 



Phard dr Tcq 

Thus the two approaches are consistent if we set 

r '' 



2Te, 



(48) 



2Te. 



(49) 



The consistency condition P^ may be further used in Eq. (24a) of Ref. [9 to 
show that in the limit |^| < 1 it agrees with our fundamental expression (IM)) . To do 
so, we rewrite Eq. (24a) in the analogous form, namely 

dx _2x^ _ ATHmsJx) 

dr T 3 

In Fig. \5\we show our function H{^) and the function Hms{0 (? = 2; — 1) multiplied 
by the factor 15/2. The two functions agree for small values of £,, hence again the 
consistency is achieved between the two approaches if the condition (|49l) is applied. 
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Since the authors of Ref. [9^ showed that their approach is reduced to the 2nd order 
viscous hydrodynamic equations of Israel and Stewart in the leading order in ^, the 
comparison showed in this Appendix shows that our approach is also reduced to viscous 
hydrodynamics in the considered case. 

We note that H{^) and (15/2)7Jms(0 differ for large values of £_. In this region 
there are no direct hints about the behavior of matter from the underlying kinetic 
theory and the form of the function i?(^) is not constrained. It may be postulated and 
eventually verified by making comparisons of the model results with the experimental 
data. Such analyses represent the most interesting applications of ADHYDRO. 
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